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Abstract 

Patterson, Schnell, &: Doyle, Nature 497 , 475 (2013) introduced a microwave experiment that 
produces a perpendicularly polarized molecular emission when chiral molecules are resonantly 
excited in the presence of a Stark Field that is then adiabatically switched off before observation of 
the emission. The sign of the signal is opposite for two stereo-enantiomers and thus the magnitude 
of the signal gives the enantiomeric excess of the sample. This paper presents a detailed presentation 
of the theory behind this and provides expressions for the absolute calculation of the expected signal 
strength for different transitions. 
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Since Pasteur recognized in 1849 that optical rotation arose from what we now call chiral¬ 
ity on a molecular scale [1] , the topic of chirality has been central to chemistry [2] . Pasteur 
also recognized that all chiral compounds uniquely arise from biological sources [3], which 
likely played a role in motivating his later work experimentally disproving spontaneous gen¬ 
eration. The relationship of molecular chirality to the underlying symmetries of physical 
laws is discussed in Barron jl] and Quack How chiral life arose from a largely achiral 

chemical soup remains one of the fundamental questions of modern science [7] . The contin¬ 
ued importance of Chirality to modern Chemistry is demonstrated by a recent search for 
’’enantiomer” in a scientihc bibliographic database which returned nearly 10,000 publications 
in the last hve years alone [8]. Despite this, the typical practicing molecular spectroscopist 
almost never directly deals with molecular chirality in his or her work because the optical 
properties of enantiomers are almost identical. The development of sensitive spectroscopic 
methods to distinguish enantiomers has long been a goal of many spectroscopists. 

In an exciting recent publication, Patterson, Schnell, & Doyle [9] presented a microwave 
spectroscopy experiment that uniquely generates signals for chiral molecules and then only 
if there is an enantiomeric excess in the sample. In their method, a chiral molecule is excited 
with linearly polarized microwave radiation resonant with an electric dipole transition while 
subjected to a Stark (DC Electric) Field oriented perpendicular to the polarization direction 
of the excitation held. In chiral or achiral molecules, this excitation produces a free induction 
decay (FID) emission held that is polarized parallel to the excitation held. However, uniquely 
for a chiral molecule, an FID is generated by the sample once the Stark field is switched off 
that is polarized perpendicular to both the Stark and excitation held directions. The signal 
is proportional, in the limit of weak excitation and Stark state mixing, to the product of 
the three projections of the dipole moment on the inertial axes of the molecule. While one 
can dehne the orientation of the inertial axes to make any two the dipole moments positive, 
the sign of the third will be hxed if one demands a right-handed coordinate system. This 
triple product of dipole moment projections has opposite sign for two enantiomers of a chiral 
molecule and is zero for any achiral molecule. For any racemic sample, the perpendicularly 
polarized FID helds from the two enantiomers will destructively interfere, again producing 
no perpendicularly polarized signal. 

More recently, Patterson et al. uniiii], Shubert at al HIE], and Lobsiger et a/. [TT] have 
reported results for a related method where the DC Stark held is replaced by a resonant 
AC held, resulting in a chiral selective four wave mixing signal. While there are obvious 
advantages of the four wave mixing approach, there are cases where the Stark mixing ap¬ 
proach may be preferred. For one, it is easier to use a microwave cavity to enhance the 
emitted power by the Q of the resonator. For another, it may prove easier to relate the 
absolute phase of the observed emission to the product of the three electric dipole moment 
components which can be used to determine which enantiomer is in excess in a sample, i.e. 
combined with predictions of Quantum chemistry, to determine the absolute conhguration 
of a chiral sample. Because the DC Stark mixing does not suffer from the same phase 
mismatching issues as the four wave mixing experiments, it may also proved easier to quan¬ 
titatively relate the observed chiral signal strengths to the size of the enantiomeric excess 
(ee) of a sample, particularly for cases where does not have a sample of the same molecules 
with known ee to use as a standard. 

In this paper, I will present an analysis of Stark-Field induced Chiral signals. [9] I will 
give results that will allow for the calculation of the size of the chiral emission, relative 
to the regular, achiral microwave free induction decay emission, for any molecule given its 
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rotational constants and dipole moment projections. We will start in the weak field limits 
where perturbation theory can be used. This will be followed by analysis of the cases of 
both strong Stark mixing and also pulses that produce maximal coherence in the sample. 


I. SIGNALS IN THE PERTURBATIVE LIMIT 

Consider a rigid chiral, asymmetric top molecule with rotation constants A, B, and C 
and projection of the molecular electric dipole along these axes of /Xq, fib and fic respectively. 
Most chiral molecules have no symmetry operations beyond the identity (point group Ci) and 
for these one expects all three components of fi to be nonzero, though obviously symmetry 
does not specify the relative sizes of the components. Some of the simplest chiral molecules, 
such as H2O2, have a C2 or higher axis of symmetry, which is necessarily along one of the 
principle axes and the permanent dipole moment is also along this symmetry axis. As such, 
the Patterson et al. [a HO] methods of chiral detection will not work for these molecules. 

The zero field rotational levels of an asymmetric top [15] can be written, in increasing 
energy order at fixed total angular momentum quantum number J, with labels | J, r, M > 
where r = and M is the projection of J on the laboratory Z axis. The more 

common Kp, Kg labels, where p is shorthand for prolate and o oblate limits, can be calculated 
using Kp = Floorand = Floor(j _ xhe field-free rigid rotor hamiltonian 
is invariant to reflection in the planes normal to the inertial axes and has symmetry of the 
Klein 4-group (Vierergruppe). The energy eigenfunctions transform under the symmetry 
operations as one of four irreducible representations that are defined by whether Kp and 
Ko are even or odd integers. The wavefunctions can be written as linear combinations of 
symmetric top functions: 

J 

\J,r,M) = AK,r\J,K,M) (1) 

K=-J 

where | J, K, M) are symmetric top functions which are functions of the Euler angles that give 
the orientation of the inertial axes relative to the laboratory fixed axes and can be expressed 
explicitly in terms of Wigner D matrices [ 12 ]. K is the quantum number for projection of 
angular momentum on the molecular fixed inertial unique axis (usually either A or C) and 
M is the projection of J on the laboratory fixed Z axis. The values of Ak,t only depend upon 
Ray’s asymmetry parameter k, = (25 —A —C')/(A — C').[T5] Likewise, the energy eigenvalues 
in zero field can be written E{J, r) = \{A -|- C)J{J -|- 1) -|- |(A — C)Ejr{K), where Ej^in) 
are the energy eigenvalues for the rigid rotor Hamiltonian with rotational constants A = 1 , 
B = K, and C = — 1 . 

Before considering the general case, it is instructive to first consider the particular case 
of driving the Oqo —t lio transition as was done experimentally by Patterson, Schnell, & 
Doyle p. We select our axis system such that the Stark Field is along the laboratory 
Z axis, the excitation pulse is polarized along the X direction and the chiral signal will be 
polarized along the Y direction. Note that this is different from what was used in [9] but this 
change has no physical consequence and I find this choice simplifies the analysis somewhat 
as then M remains a good quantum number. We neglect AJ = ±1 mixing terms in the 
Stark effect. With this approximation, there are no Stark contributions to the Oqo state. 
The lio state has Stark Matrix elements with the In state with an a polarized dipole and 
energy spacing B — C and state Iqi with a h polarized dipole and spacing A — C. We will use 
the notation |J, r, M)^ to denote the energy eigenstate in a Stark Field that adiabatically 
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connects with the held free |J, r, M) state. To hrst order in the Stark Field we have: 


|1io,m)^ — |1io,m) — 


B-C 




A-C 


|1o1,m) (1o 1,M|0bz|llO,M) (2) 


where are the direction cosine matrix elements between the molecular hxed axes a = 
a, b, c and the laboratory hxed axes G = X,Y, Z. The sample is driven by an excitation 
pulse, EexCos{ut), that is on resonance with the degenerate |Ooo;o)^ —t |lio,±i)^ transitions. 
For times after the end of a weak excitation pulse of length At, and including only near 
resonance contributions, the wavefunction is 

|^(t) = |Ooo)'+ ^ om |li0;M)'exp (-iE'(lio)t/h)) 


M=±l 


aM=±i — 


iEe^At 

2h 


iEe,^At 

2h 

iEe^^At 


he (llO,Ar|0cX:|Ooo)-(llO,M|0az|lll,M) (lll,M|0bx|Ooo) 


Est^at^b 
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B-C 


(1iO,m|06z|1o1,m) (loi,M|0ax|Ooo) 


l^c 






A.2i^c 


B-C 
Est^at^bi,A B'j 
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i — 
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(3) 


4v^h L {A-C){B-C)_ 

Note that the contributions from mixing of the In and Iqi states with the lio state are 
opposite in sign. We assume that the Stark Field is turned oh adiabatically before the FID 
of the sample dephases. For detection with held oh, we have polarization in the Y direction 
of 


By = 


AN ^ a(M) exp {-iE{lio)t/h) (Ooo|hc0cy |1io;m) + c.c. 


M=±l 






I AN 


EsEe^At 


±2inc - 

A-B 


12h {A-C){B-C) 


Est^at^bi^A B^ 

{A-C){B-C) 


he - 


1 

V6 


+ c.c. 


(4) 


where AN = Y(Ooo) —-Y(lio)/3 is the thermal diherence in number density for the transition 
and c.c. stands for complex conjugate of what comes before. As predicted, the polarization 
to zero order in E^ is zero because of the cancelation in the sum over M = ±1. The 
Stark induced signal is proportional to the product of hahfehc and thus is of opposite sign 
for two enantiomers. Thus, if we have a mixture of enantiomers, the net amplitude of the 
polarization will be proportional to the diherence in the number density of the R and S 
forms, i.e. proportional to the enantiomeric excess. We can compare this signal to the X 
polarized polarization that leads to the standard FID; 


Px = 


AN ^ a(M)exp(-iF;(lio)t/h) {Ooo|hc0cx|lio;M)+C.C. 


M=±l 




M=±l 


4v^h 


±2ific 


{A-C){B-C)\ 


he T 


Ve 


+ c.c. 


6h 


( 5 ) 


4 



Both Py and Px are phase shifted by 90 degrees from the excitation wave. The FID held 
emitted these polarizations will be further shifted by 90 degrees, so, both net helds will 
be out of phase with the driving held, should it had been continued into the detection 
interval. This result can be physically understood if one considers an experiment that seeks 
to simultaneously detect both Px and Py. Imagine that after the sample (in the direction of 
propagation of the excitation wave, which is Z) we place a horn capture the X polarized FID, 
and then we rotate it to detect the Y polarized wave. The sign of this wave (relative) to the 
X polarized wave will change depending upon if we do a clockwise rotation (looking in the 
direct of the wave propagation) or counter-clockwise rotation. In the clockwise rotation case, 
we are dehning Y to form a right-handed coordinate system with respect to the previously 
dehned X and Z axes, while the opposite rotation dehnes a left-handed coordinate system. 
The net FID has rotated like in normal optical rotation. 

Before we turn off E^, we have to calculate the transition of |Ooo) to the Stark Field mixed 
state. That leads to additional contributions to Py: 


Py 


AN 

M=±l 

-^Shah'fe 

~ A-C 

AN 

M=±l 

Esf^af^b 

~ A-C 


fj'c (Ooo|0^ 


cVl 


)- 


B-C 


(Ooo| 0 bv|lll,Ar) (lll,Ar|0az|llO;M) 


(Ooo | 0 ay I 1o1,m) (1o1,M 106 Z I 1i0;m) 

f 1 ^ EsHa^b ( i \ 

V Ve; B-C \ ^) 



+ c.c. 
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substituting in from above for a{M) gives 


Py = 


iAN EexAt 


iAN EexAt 




E(=F1) 


M=±l 

-l-c.c. 

M=±l 


P-c T ^ 


Pc + 


. EsHafJ^biA - B) 


2{A-C){B-C)\ Y 


EsfJjafJ^bi^A 5 ) 

2{A-C){B-C) 


Hc±i 


■ Ejjfjjdfibi^A f?) 


2{A-C){B-C)\ 

(7) 


-I- c.c. = 0 


and thus demonstrating that the Y polarized signal (i.e. perpendicular to the excitation 
polarization) does not appear until the Stark held is changed. Note that if we had inverted 
the direction of the Stark held, which changes the sign Es, we would double the amplitude 
of the chiral FID compared to simply shorting the held to zero. 

Before considering the ehect of excitation on a generic rotational transition in the presence 
of a Stark Field, let us consider the held free excitation of a transition |J",r") —)■ 
with X polarized radiation. We assume that we start with a thermal population with states 
IJ, r, M) equally populated for hxed J, r, with total population in a level J, r of N{J, r). The 
excitation pulse will be assumed to have electric held amplitude E^^, duration At and an 
angular frequency detuning from exact resonance Au. a = a,b, or c are the directions of the 
transition moments in the principle axis system of the rigid molecule. Following excitation. 
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we will generate a sample polarization in the Y directions 


Py = AN 




^—iAujAt 


) 


2hAu 




(J", r", M'Va.rl^', r', M') (J', r', M'|0a,x| J", r", M") + c.c. (8) 

a,M",Ar'=Ar"±l 


where ojq = {E{J', t') — E{J”, is the angular frequency of the emission and 


V 2J" + 1 2J' + 1 ) 


(9) 


is again the thermal population difference between states of hxed M. We can decompose 
the direction cosine matrix elements as: 


= (J, r|0„| J', r') (J, M') (10) 

(J,r|0„|J',r') = (J|0j|J') AK,rAK'y {J, K\<P^\J', K^) (11) 

K,K' 


where the matrix elements of 0j, and 0 j,m are dehned in Table 4.2 of Townes & 
Schawlow na. For a given transition, at most one value of a will give nonzero matrix 
elements M''\(j)a^G\J' n', M'). We can use the symmetries 


(J', r', M"±l y I 


(J', t\ M" ± l\ct>^,x\J'\ r", M") (12) 

± l|0„,x|J",r",M") (13) 


to show that 


(J", r", r', M" ± 1) (J', r', M" ± l|0„,v| J", r", M") = 

- (J", r", I J', r', -M" ^ 1) (J', r', -M" ^ . r", -M") (14) 

Thus, as long as the population of the states with M" = +M and M" = —M are equal, the 
terms with M" = M,M' = M ±1 will cancel those with M" = -M, M' = -M =f 1 and we 
get no net polarization perpendicular to the excitation polarization and thus no Y polarized 
contribution to the Free Induction Decay (FID) emission. This result is independent of the 
labels we give to the axes of excitation and detection. Physically, this is a consequence of the 
fact that the sample is not birefringent if we have zero orientation of the sample. Even achiral 
molecules in states of dehnite M are optically active, i.e. have different susceptibilities for 
left and right handed light. Such individual states have what Barron calls “false chirality” |1] 
but not chirality, based upon the symmetry properties under time reversal. 

Consider application of a DC Stark held along the Z axis with electric held Es, and for 
now assume it is small enough that we can use perturbation theory. To hrst order, we get 
for the wavefunctions in the presence of the held: 

|J,r,M >' 


= I J, r, M > 

. 7+1 

- E EE 

Js=J-l Ts a 


(15) 




{E{J,t)-E{J,,Ts)) 


\Js,Ts,M) 
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We will now consider excitation on a transition \ J'' ,t'' , M")' —)■ \J' ,t' , M')' with, as 
before, a resonant field with amplitude polarized along the X axis that is nonzero for 
time interval At and that the excitation is sufficiently weak we can use hrst order time 
dependent perturbation theory. Following the excitation pulse, the wavefunction for each 
initially Stark perturbed rotational state \ J" ,t" ,M")' can be written: 






+ a{M”, M') IJ', r', M')' 

M' 


(16) 


Making the rotating wave approximation, we can write for times after the end of the exci¬ 
tation pulse that: 




2h{Au{M',M")) 


/3 


(17) 


where Au{M', M") is the M dependent detuning for the transition. Due to the perpendicular 
polarization, we have nonzero contributions only for M' = M" ± 1. To simplify things, we 
assume that the second order Stark splitting of the transition is sufficiently small that 
Au{M', M")At << 1, in which case we can treat all M components as on resonance. This 
allows us to write 

rp 

a(Ar, M') = ^ {J', r', Af'l' .^^.v \J", t", M")' (18) 


More generally, we can replace E^xAt by / Eex{t)dt. We will keep only terms up to linear 
in Es. This allows us to write for the transition matrix elements 


^ (f, t', Af'l' \j", t", M")' = 5^ {J', t', Af'l 4>f.x \r, t", M") 

P /3 

' (J', r', M'l I M') ( M '\| J", r", M") 


— Es 

a,p,Js,Ts 


E{J',t')-E{Js,Ts) 

(J', r', M'l I J., D, M") (M"| | J", r", M") 

E{J",t")-E{Js,Ts) 


(19) 


Since we adiabatically turn off the Stark Field before detecting the emission polarized along 
T, we calculate the emission with |J, r, M)^ replaced by |J, r, M), i.e. the zero order wave- 
functions. This leads to a molecular polarization polarized along T, Py, that results in 
emitted FID held oscillating at the excitation frequency and polarized along Y, 


Py = AN Y + C.C. (20) 

Above, it was demonstrated that the terms zero order in Eg give zero contribution after 
the sum over M levels. The terms linear in Es give: 

—iEsEexAtAN 


Py = 


2h 


/i«/i/3h7 {J''x"\(l)y\J',T') (j",M"|0y|j',M') X 


a,/3,7 M',M",Js,Ts 

(J',r'|0^|J„r,) (J',M'|0z|J.,M') (J", r") (J„ M'|0x| J"M") 

P(J'r')-P(J„r,) 

(J'r'|0;3|J.,r,) (J'M'|0x|J.,M") (J„ r,|0„| J"r") (J„ M"|0z| J", M") 


P(J",r")-P(J„r,) 


-I- c.c. 


( 21 ) 
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Given the selection rules for an asymmetric top, we can have a cycle of transitions (J", r") —)■ 
{Js,Ts) —)■ — )■ only if a 7 ^ 7 ^ 7 , in which case a,/3,'j are a permutation 

of a, b, c. The Py signal will only be nonzero if all fia, l^b, he 7^ 0, which is true only for 
molecules of Ci symmetry, which are chiral. Only M' = M" ± 1 give nonzero contributions. 

We can do the sums over M quantum numbers if we specify the J quantum numbers. 
For an R branch transition J" = J, J' = J + 1, J* = J, J + 1. we get after summing over M 
quantum numbers: 


Py{J,t'' ^ J + l,r') = 

2EsE,,AtAN ^ + 1)(2J + 3) ^ (Jr"|0^| J + 1, r) 

Oil 

J+2 J+2 


- (J + l,r'|0a|J + l,Ts) {J + l,Ts\(l)f}\J,T”) ^ 
+ (J+ l,T'\4>a\J,rs) {J,Ts\4>/i\J,T'') 


J 


Ej+iy — Ejyi^r, 

J 


F/jj-i t' — E, 


J,Ts 


E T -r" — E 




Ejy — Ejyi^Ts 
+ c.c. (22) 


For a Q branch transition with J' = J" = J, Jg = J — 1, J, J + 1. Summing over the M 
quantum numbers gives the following result 


Py{J,t" ^J,t') = 
2EsEe^AtAN (naUb/^c) 
3h 


e-'““V(J+l)(2J + l) ^ (Jr"|0^|J,r') 

a,l 3 ,'y,Ta 


X 


{J,T'\(j)a\J + l,Ts) {J + l,Ts\(j)fi\J,T'') 



2J + 3 

J,t' — Ejyi^rs 


2J + 3 


E T -r" — E 


J+ljTs 


- (J,r'|0a|^,r,) (J,r,|0/3|J,r") (- - 

- (J,r'|0a|7 -It's) {J - 1,Ts\4>^\J,t”) 


1 - ^\ 

Ejys Ej y" Ej^^^ j 


f 2J-1 

2J-1 Y 

1 

1 

J) 

1 

Ej^t" j _ 


+ c.c(23) 


The triple product of direction cosine matrix elements will be imaginary, which means 
that the FID in the Y polarization will be phase shifted by ±7r from the driving radiation, 
like the X polarized polarization, so both will lead an emitted held out of phase with the 
driving radiation. Each of the direction cosine elements (^q) scales ~ J“^, so except for the 
J" = J' = Jg terms, we get a net ~ J scaling of the terms, which rehects the degeneracy 
of states. In most cases, the dominant terms will be Stark Mixing of asymmetry doubles as 
these have the smallest energy differences. For such mixing, the chiral signal on Q branch 
transitions will be weak relative to R branch transitions. 

It is worth noting that the sum over M quantum numbers in the triple product of direction 
cosine factors gives zero unless the three laboratory projects are a permutation of X, Y, Z, i.e. 
that the polarizations of the Stark, Excitations, and Detection helds are each perpendicular. 
This reflects the fact that the chiral signal arises from projections of the three dipole moment 
components from the rotating molecular frame onto a laboratory frame dehned by the three 
electric held directions. If we project the three vectors onto a plane, then there is no way 
to determine if we are projecting from a left-handed or right-handed axis system formed by 
the three dipole projections on the molecular axis system. 



While the Stark Field is unchanged from its value during the excitation, the Y polarization 
is calculated using the Stark perturbed wavefunctions: 

Py = AN + C.C. (24) 

This leads to additional terms which are almost the same as above but when summing over 
M quantum numbers we get terms 

(J", M"|0y IJ', M" ± 1) (J', M" ± l\cj)x\Js. M") ( M" ± l|0z| j", M") + 

M" ,M'=M"±1 

± 1) ± l|0yI J„M") (J„M" ± l|0z| (25) 

and 

(J", M"|0y 1^., M" ± 1) ( M" ± l|(^z| J', M" ± 1) (J', M" ± l|0x| J", M") + 
± 1) (J„M" ± J',M" ± 1) ± l|0yI (26) 


where in both cases, only the hrst term in each sum is present when we detect with the 
Stark held off. Using Eq. 14, it is evident that the sum of the hrst and second terms in the 
above equations are identically zero, i.e. no signal exists until the Stark Field is changed and 
the amplitude of the FID will be proportional to the change in Eg. We could also do the 
excitation with Eg = 0 and the detection with hnite Eg. This will produce a signal of the 
same magnitude but changed in phase. The signal could be doubled in size by changing the 
sign of Eg between excitation and detection. However, in these cases, the FID will contain 
dephasing from the diherent M components of the transition. In many cases, sufficient 
microwave power is available so that the excitation time At will be signihcantly shorter 
than the dephasing time of the field free FID, and in these cases it is better to have the 
Stark dephasing with the shorter excitation time instead of the longer detection time. 

If we restrict the Stark mixing to only AJ = 0 coupling terms, the signal (up to scale 
factors) depends only on the asymmetry factor k. Numerical calculations demonstrate that 
this chiral signal is nonzero for all R branch transitions of the chiral molecule for typical 
values of k. As could be expected, the predicted signal is particularly larger for transitions 
to levels with small asymmetry splittings as then the Stark mixing is larger due to the small 
energy denominator. In these cases, the pairs of transitions to the two states that make up 
the asymmetry doublet have nearly equal and opposite sign signals. 

Because the enantiomer identihcation depends upon the phase of the Y polarized FID 
signal, one should account for the fact that the Stark shifts will lead to some dephasing and 
phase shifts of the FID that depend upon the time delay and period for Eg to be ramped 
down. One also needs for the time dependence of Eg to be sufficiently slow for the states to 
evolve adiabatically. This likely puts some limit on the magnitude of the Stark mixing. If 
the Stark interaction is such that one gets substantial state mixing, then the Stark shifts are 
comparable to the energy levels splittings and the two conditions of small dephasing of the 
M components and that one has adiabatic evolution of the rotational states will be incom¬ 
patible. Table I contains the predicted chiral signals for the R branch transitions calculated 
for the case that k = —0.706 which is the value appropriate for 1,2-propanediol ini. the 
molecule used by Patterson, Schnell & Doyle to experimentally demonstrate this effect [U]. 
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The column labeled P{Y) should be multiplied by (^ ^ 2 j"+i ^ ~ ^j 

to put the numbers in absolute units. The column labeled S has the line strength factor 
for each transition, S = w a I {J'y'’''^ , M”) p and is one of the factors 

contributing the strength of the X polarized FID emitted by the sample. 

II. HIGH STARK FIELD MIXING 

The energy spectrum of an asymmetric top consists mostly of nearly degenerate levels 
known as asymmetry doublets. These correspond semiclassically to symmetric and anti¬ 
symmetric combination of states rotating in opposite directions around either the a or c 
axes and the splittings of asymmetry doublets can be well approximated by semiclassical 
tunneling integrals between the semiclassical trajectories corresponding to classical rotation 
around these axes [H] . For such states, the Stark effect at low held is dominated by mixing 
of these asymmetry doublet states if the corresponding dipole component is nonzero, which 
is the case of the chiral molecules we are considering in this work. We will now assume we 
have an a type asymmetry doublet and we will select the direction of the a axis to make 
Pa > 0 which we are free to do. The results apply to c type asymmetry doublets if we invert 
the a and c labels. 

We consider the transitions \ J,±K) \J + l,±(iF -|- 1)). We will use symmetric top 
functions. We write the asymmetry splitting in the \J,±K) state as 2Ai/" and that in the 
\J -|- 1, ±(iF -|- 1)) as 2Ai/'. Let uq = {B + C)J (2A — B — C){K 1/2) be the transition 
frequency without asymmetry splittings. In zero Stark held, we have four transitions at 
uo - XE"/h - AE'/h {Jk,j-k J + ^K+i,j-K+i-, b type), z/q — AE"/h + AE'/h {Jk,j-k —t 
J + 1k+i,j-k,c type), vq AE"/h - AE'/h {Jk,j-k+i J + Ik+i,j-k+iX type), and 
Uo -h AE"/h + AE'/h {Jk,j-k+i J + lK+i,j-K,b type). 

As we turn up a Stark Field, each of the two asymmetry doublets will split from each 
other. Initially, we will get second order energy shifts of ±[fiaKMEs/J{J + l)]^(2Ai?)“^, 
with -|- for the upper asymmetry component and — for the lower. As the Stark shifts grows 
relative to the zero held asymmetry splittings, the shifts in each doublet evolve to hrst order 
in Es with energy values hujKM = hujKTf^aKMEs/J{J+l), where hujK is the mean energy 
of the J, K asymmetry doublet. The transitions that adiabatically evolve from the b type 
transitions will loose intensity, going to zero intensity at sufficiently high Es- The transitions 
that evolve from the c type transitions in zero held will strengthen, taking up the decreasing 
strength of the b types. In high Stark Field (i.e. when \fiaEsKM\/J{J + 1) >> AE', AE"), 
we can ignore the asymmetry slitting and the eigenstates will be nearly pure symmetric 
top states as the diherence in Stark energy will quench the tunneling. Consider excitation 
on the Stark component \J,K,M) —)■ | J-|-1, iF-|-1, M-|-1) {K,M > 0). This will be 
exactly degenerate with the transition \ J, —K, —M) —)■ | J -|- 1, —K — 1, —M — 1). For high 
|M|, these are the strongest transitions in the Stark multiplet of the J,K transition. If we 
polarize these transition with X polarized radiation, we will create excited state amplitudes 
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proportional to: 


a{±K,±M) = 


iEe^At 

2h 


(< J + 1, + 1, M + + f^c(l^c,x\'J, K, M >) (27) 


bjK = 


7 r r -1 

- ^ - OjkOjm [IJ^b =F ll^c\ 

(J + 77 + l)(J + 77 + 2)~ 

2(^J + 1) (2J +1)(2J^ “l“3) 


(28) 


The states that \ J, AK, AM) and | J + 1, A{K + 1), ±(M + 1)) will adiabatically evolve 
to, as the Stark Field is reduced to zero amplitude, depends upon the sign of KM. 
We will assume here that KM > 0 in both initial and final quantum states (which 
are the stronger transitions), which implies that the states will evolve to the lower 
asymmetry states, which are K^AM) — \ J — K^AM)) for the lower state and 

(I J + 1, 77 + 1, ±(M + 1)) — I J + 1, —K — 1, ±(M + 1))) for the upper state. The lin¬ 
ear combination of the two lower states will give for X and Y polarized emission 


Px = -iJ,c{ihjKK){AbjM)a{,AK, AM)e + c.c. 

± 

= ,4 + C.C. (29) 

Py = -/ic(i&ji^/ie)(^7jM)a(±^, ±M)e~'‘^°* c.c. 

± 

= + c,c (30) 


The strongest transitions 
factor b'jj^b'j^ becomes 


for perpendicular excitation are when J 



2J+ 1 

4(2J + 3) 


K 


M for which the 
(31) 


We thus see that the amplitude of the Y polarized emission saturates at high Stark 
held such that the ratio of the amplitudes of Y polarized emission to X polarize emission is 
proportional to fib/K- If had taken KM < 0, the pumped states will adiabatically evolve 
to the upper levels of the asymmetry doublet as the Stark Field is turned off. This will give 
a zero held Y polarized molecular emission with the same expression as above except for a 
sign change. Note that since we have selected the coordinate system to make fia > 0, this 
phase of the Py emission will remain opposite in sign for a pair of enantiomers and thus its 
size will proportional to the e.e. of the sample. 

There will likely be a substantial phase shift during the time that E^ is turned oh, but 
this should ehect Px equally as Py and thus the relative phase of the these signals will give 
the chirality. Figure 1 shows a plot of the relative excitation amplitudes for the b and c type 
transitions of J” = K" = M" = 2 —)• J' = 77' = M' = 3 transitions of propanediol, which 
has a center transition frequency of 46 GHz as a function of 77s, using the constants given 
by Lovas et. al. m 
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Transition Intensity vs. Stark Field 



Stark Field in V/cm 


FIG. 1: This shows the intensities of the transitions between a pair of asymmetry doublets as a 
function of Stark Field. Parameters used for J” = K'^ = M" = 2 —)■ J' = K'^ = M' = 3 transitions 
of propanediol 


III. EFFECT OF VARIATION IN EXCITATION STRENGTH 


The above results were derived in the limit of weak excitation, we will now consider 
polarization pulses that move significant population between the two levels of a transition. 
If Eq is sufficiently large that the individual Stark components M” —)• M' are resolved 
in the excitation, then both Px and Py FID signals will be proportional to sin(0) where 
6 = {J', K', M'\'^a4>a,x\J'', K", M")' /2h is the {Eg dependent) excitation flip 

angle. Thus, FlD’s of both polarizations will be optimized for x 12 excitation. In the high 
Stark field limit considered above, we will have a flip angle equals |a| as given by Eq. 27 
Note, however, that this is for a single Stark component. 

If Es is sufficiently weak that the transitions \J" ,t" , M") —)■ |J',r',M"±l) are both 
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excited, then the pulse creates a coherent superposition of the |J',r',M"±l) states. 
We hrst consider the Stark Field free case where we are pumping \J",K",M'') —)■ 
\J', K', M' = M" ± 1) transitions. Following an X polarized excitation pulse we have 


= cos ^ ^ -^- +*sin ( —— 1 e-^— x 


(32) 



|J',r',M" + l) (J',M" + l|0x 

\r,M") + \. 



VI + l|0x| 

J", M") 

P + 

\{J',M"-l\(j)x\ 


p y 


with M” dependent Rabi frequency given by: 

V a 

Vl (J', M" + l\(t)x\J", M") |2 + I (J', M" - l|0x| J", M”) |2 (33) 

We can write 


+ l\(t)x\J\M") |2 + I - l\(t)x\J\M") |2 = 

2(J" + 1)(J" + 2) + 2M^ if J' = J” + I 

2 J"( J" + 1) - 2M2 if J' = J” (34) 


In the limit that we can treat the Stark Mixing with hrst order perturbation theory, we 
can write the Rabi frequencies linear in Stark Field as 


u{M'\Es) = u%M”) - Es5u{M”) 


(35) 


6uj{M") = (36) 

^ ^exhah/3 (^^ M" + l\4>x\J", M") 6+,.,/?(M") + ( M" - l|0x| J", M") &-,.,/?(M") 
f- n VI (J', M" + l|0x| J", M") p + I (J', M" - l|0x| J", M") |2 


6±.a,e(M'') = 5^ 

Js-,Es 


(J', r', M" ± 1| I M" ± 1) (M" ± 1| ^^^x | J", r", M") 

E(J',r')-^(J.,r,) 

(J', r', M" ± 1| Vs, M") (J„ M"| | J", r", M") 


^ E(J",r")-^(Js,r,) 

The polarization produced parallel to R will be: 


(37) 


Py = 'E\ (J", r", x (38) 

M",-f 



|J',r',M" + l)(J',M" + l|0x 

|J",M") + |, 



VI V',M" + l|0x| 

J", M") 

|2 + 

|(J',M"-l|0x| 


p y 


+ C.C. 
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To first order in Es 

sin(a;(M", Es)At) = sin — Es cos 6uj{M'')At (39) 

As before, the term zero order in Es generates no Y polarized FID as long as the population 
of states with M" = +M and —M are equal. 

^ cos(a;°(M")At) 

nj', t', M" + 1 ) (J', M" + l\(j)x\J", M") + IJ', r', M" - 1 ) (J', M" - l\(j)x\J\ M") 

X ((/, M" + l|0x| J", M") 6 +,„,/ 3 (M") + (J', M" - l|0x| J", M") 6 -,«,/ 3 (M")) + c.c. (40) 

IV. CONCLUSIONS 

The breakthrough experiment of Patterson, Schnell & Doyle has demonstrated a 
novel way to study chiral molecules using microwave spectroscopy, which already allows 
(through isotopic substitution) the highest resolution way to determine equilibrium structure 
of isolated molecules but previously was not able to determine the chirality of a molecule. 
This paper presents a thorough analysis of the theory behind this experiment and presents 
expressions that can be used to calculate the experimental observables in this and related 
experiments. 
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TABLE I: Calculated Chiral FID for k = —0.706 


Kp ". Ko " 


J Kp '. Ko ' 

P(Y) 

S 

polarization 

O 

o 

o 


1 0,1 

-0.014345093 

0.333333333 

a type 

o 

o 

o 


11.1 

0.6651107 

0.333333333 

b type 

o 

o 

o 


1 1.0 

-0.48409894 

0.333333333 

c type 

1 0,1 


2 0,2 

-0.093857917 

0.663570114 

a type 

11,1 


2 0,2 

-0.183313521 

0.207471598 

b type 

1 1.0 


2 0,2 

0.129373065 

0.128958289 

c type 

1 0,1 


2 1,2 

0.248745836 

0.5 

b type 

11,1 


2 1,2 

0.176560314 

0.5 

a type 

1 0,1 


2 1.1 

-0.103834802 

0.5 

c type 

1 1.0 


2 1,1 

-0.104109539 

0.5 

a type 

1 1.1 


2 2,1 

7.322764362 

0.5 

c type 

1 1,0 


2 2,1 

7.042634499 

0.5 

b type 

11.1 


2 2,0 

-6.442956114 

0.459195069 

b type 

1 1.0 


2 2,0 

-7.832401565 

0.537708378 

c type 

2 0,2 


3 0,3 

-0.312002831 

0.988376378 

a type 

2 1,2 


3 0,3 

-0.170462573 

0.464492753 

b type 

2 1,1 


3 0,3 

0.057202463 

0.211955453 

c type 

to 

o 


3 1.3 

0.230292826 

0.69283623 

b type 

2 1,2 


3 1,3 

0.425185158 

0.887727716 

a type 

2 0,2 


3 1,2 

-0.05625524 

0.640647363 

c type 

2 1.1 


3 1,2 

-0.135165681 

0.887417203 

a type 

2 1,2 


3 2,2 

1.051541456 

0.555555556 

c type 

2 1,1 


3 2,2 

1.104838224 

0.555555556 

b type 

2 2,1 


3 2,2 

0.010202148 

0.555555556 

a type 

2 1,2 


3 2,1 

-0.841909744 

0.424396136 

b type 

2 1.1 


3 2,1 

-1.249563689 

0.676933435 

c type 

2 2,0 


3 2,1 

-0.015705437 

0.555777842 

a type 

2 2,1 


3 3,1 

168.1415563 

0.847860635 

c type 

2 2,0 


3 3,1 

152.5967058 

0.816611135 

b type 

2 2,1 


3 3,0 

-152.1885753 

0.81455008 

b type 

2 2,0 


3 3,0 

-168.4892139 

0.849764516 

c type 
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TABLE II: Calculated Chiral FID for k, = —0.706, continued 


J” Kv " . Ko " J’ Kv '. Ko ' 

P(Y) 

S 

polarization 


3 0,3 ^ 4 0,4 

-0.783932511 

1.30828008 

a type 


3 1,3 ^ 4 0,4 

-0.286248405 

0.772842846 

b type 


3 1,2 —^ 4 0,4 

0.035655368 

0.245690525 

c type 


3 0,3 ^ 4 14 

0.343979911 

0.926718354 

b type 


3 1,3 ^ 4 1,4 

0.892992397 

1.246446683 

a type 


3 0,3 ^ 4 1,3 

-0.047856102 

0.738091198 

c type 


3 1,2 ^ 4 1,3 

-0.125335263 

1.24461383 

a type 


3 2,2 ^ 4 1,3 

-0.257571118 

0.187451354 

b type 


3 2,1 ^ 4 1,3 

0.247715437 

0.162755388 

c type 


3 1,3 —^ 4 2,3 

0.278880711 

0.613967307 

c type 


3 1,2 ^ 4 2,3 

0.352860332 

0.638828886 

b type 


3 2,2 ^ 4 2,3 

0.035371194 

0.998783251 

a type 


3 1,3 ^ 4 2,2 

-0.214634099 

0.3637627 

b type 


3 1,2 ^ 4 2,2 

-0.367408049 

0.866105215 

c type 


3 2,1 ^ 4 2,2 

-0.048336179 

1.000693648 

a type 


3 2,2 ^ 4 3,2 

19.19218407 

0.918062199 

c type 


3 2,1 ^ 4 3,2 

16.89320368 

0.823592697 

b type 


3 3,1 ^ 4 3,2 

-0.000113317 

0.584011457 

a type 


3 2,2 ^ 4 3,1 

-16.65395456 

0.812548646 

b type 


3 2,1 ^ 4 3,1 

-19.38320717 

0.928278022 

c type 


3 3,0 ^ 4 3,1 

-0.001708194 

0.584015534 

a type 


3 3,1 ^ 4 4,1 

3806.704117 

1.181293203 

c type 


3 3,0 4 4,1 

3469.978829 

1.150427611 

b type 


3 3,1 4 4,0 

-3469.674136 

1.150332728 

b type 


3 3,0 —^ 4 4,0 

-3806.963474 

1.181380843 

c type 
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